Spin nematic correlations in bilinear-biquadratic 5 = 1 spin chains 
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We present an extensive numerical study of spin quadrupolar correlations in single and coupled 
bilinear-biquadratic spin one chains, using several methods such as Exact Diagonalization, Density 
Matrix Renormalization Group and strong coupling series expansions. For the single chain we clarify 
the dominant correlation function in the enigmatic gapless period-three phase for 9 £ (7r/4, tt/2), 
which is of spin quadrupolar nature with a period three spatial structure. Then we revisit the 
open problem of the possible existence of a ferroquadrupolar phase between the dimerized and the 
ferromagnetic phases. Although an extended critical region is in principle compatible with the 
numerical results, a scenario with a huge crossover scale is more plausible. Finally we study the fate 
of the dimerized phase upon coupling two chains in a ladder geometry. The dimerized phase rapidly 
vanishes and an extended gapped phase takes over. This gapped phase presumably has dominant 
short-ranged ferroquadrupolar correlations for 9 £ (— 3-7T, 4, — tt/2) and - suprisingly - seems to be 
adiabatically connected to the plaquette single solid phase of the Heisenberg S — 1 ladder and 
therefore also with the Haldane phase of isolated chains. 

PACS numbers: 75.10.Jm, 75.10.Pq, 75.40.Mg, 75.40.Cx 



I. INTRODUCTION 

The recent experimental demonstration^ of the transi- 
tion from a superffuid state to a Mott insulating state of 
atoms in an optical lattice has opened the way to novel re- 
alizations of effective quantum lattice models with widely 
tunable control parameters. Quantum magnetic systems 
can be realized by spinor atoms in an optical lattice, e.g. 
23 Na with a total 5=1 moment. Confining 5=1 
atoms to an optical lattice there are two scattering chan- 
nels for identical atoms with total spin 5 = 0,2 which 
can be mapped to an effective bilinear and biquadratic 
spin interactior 
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where we adopt the standard parametrization Jbi = 
cos 9 and Jb q — sin 9. In one dimension, the bilinear- 
biquadratic spin-one model has a rich phase diagram (see 
Fig. ^) with some well established phases: the Haldane 
gap phased, a dimerized phased, a ferromagnetic phase 
and some less well understood phases: a critical phase 
with period three correlations^, which we will charac- 
terize as having dominant spin nematic correlations, and 
possibly a gapped spin nematic phase between the dimer- 
ized and ferromagnetic phase 8 , which however remains 
controversial. 

On the square and the simple cubic lattice, the 
bilinear-biquadratic spin-one model is well understood 
for the case Jb q < 0. It exhibits a ferroquadrupolar spin 
nematic phase for 9 E (— 37r/4, — 7r/2)^i2iii. Adjacent to 
it are an antiferromagnetic Neel phase and a ferromag- 
netic phas e 9 : 10 . The region of purely antiferromagnetic 
couplings Ju , Jbq > remains to be understood, and 



could possibly contain a spin liquid region^. Recently 
the Hamiltonian JTJ on the triangular lattice attracted 
some intereslji 3 ^*^, as a possible explanation of the un- 
conventional magnetism of NiGa2S4 16 . 

For a single chain Chubukov 8 suggested the existence 
of a gapped, nondimerized phase showing dominant spin 
nematic correlations close to the ferromagnetic region of 
the phase diagram. Subsequent numerical workii could 
however not substantiate this claim and it was there- 
fore believed for a while that the dimerized phase would 
extend up to the ferromagnetic phase boundary. Recent 
quantum Monte Carlo calculations 18 and field theoretical 
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FIG. 1: Phase diagram of the bilinear-biquadratic spin-one 
chain. The firmly established phases are the Haldane, the 
ferromagnetic and the dimerized phase. We characterize the 
extended gapless phase 7r/4 < 9 < tt/2 by having dominant 
k = ±27r/3 spin quadrupolar correlations. The possible oc- 
curence of a spin nematic like phase close to — 37r/4 is inves- 
tigated and critically discussed. 
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work 19 suggest that this picture might need to be recon- 
sidered, especially in the light of possible experimental 
verifications in Bose-Einstein condensate systems^. In 
the meantime there has been a considerable number of 
publications in favor or against a spin nematic phase close 
to 9 = -37r/ l 2 ?i 21 i 22 i 2 ?i 24 , but leaving the final answer 
still open. 

The aim of the present paper is twofold. First we 
give a characterization in terms of the dominant corre- 
lation function for the extended gapless critical phase, 
whose mere existence is well established. We find in 
this case that the predominant correlations are of spin 
quadrupolar nature, with a wavevector of ±27r/3. The 
second aim is to shed some more light on the region 
9 6 [— 37r/4, — 7r/2] of the single chain. We present ex- 
tensive numerical simulation results based on different 
numerical methods and discuss possible interpretations. 
Due to the difficulty in obtaining a firm conclusion re- 
garding the spin nematic phase, we have also studied 
bilinear-biquadratic ladders. In this case the dimerized 
phase plays only a minor role in the phase diagram and 
we report an extended gapped spin nematic phase, very 
close in spirit to the proposal by Chubuko\«& for the single 
chain. Surprisingly this phase seems to be adiabatically 
connected to the well known Haldane phase of the iso- 
lated chain. 

The outline of the paper is the following: in sec- 
tion [H] we show that the established gapless phase 9 e 
(+7r/4, 7r/2) has dominating spin quadrupolar correla- 
tions, a characterization which was lacking before. We 
then move on to the region 9 € [— 37r/4, — 7r/2] and dis- 
cuss the possible existence of a spin nematic intermedi- 
ate phase between the dimerized and the ferromagnetic 
phase. We report several anomalous physical proper- 
ties encountered upon approaching 9 — > — 37r/4, which 
could be interpreted as a phase transition to a spin ne- 
matic state. However, due to numerical limitations it 
has remained elusive to pinpoint such a phase transi- 
tion, and an alternative scenario where a a very large 
crossover scale emerges as one approaches the SU(3) 
point at — 37r/4 becomes more plausible. In section 17x71 
we introduce an extension of the single chain model by 
coupling two chains to form a ladder. We show that the 
dimerized phase gives very rapidly way to an extended 
gapped phase. This short-ranged ordered phase encom- 
passes the Haldane phase at 9 — as well as its ladder 
extension, the so called "Plaquette Singlet Solid" stated, 
and crosses over to a gapped spin nematic state close to 
9 -> -3?r/4. 



II. THE ENIGMATIC PERIOD 3 PHASE 

The existence of an extended critical phase in the in- 
terval 9 € [+7r/4, it/ 2) has first been discussed in Ref. 0. 
Later work in Refs. |2ffl27>2Sl2<l3rj agreed on the gapless 
nature of the phase due to soft modes at k = 0, ±27r/3. 
The special point 9 = +7r/4 has an enlarged SU(3) sym- 



metry and is solvable by Bethe ansata^i^SiS (Uimin-Lai- 
Sutherland model), proving the existence of soft modes 
at k = 0, ±27r/3 as a rigourous result. 

Despite the consensus on the presence of an extended 
gapless phase, the nature of the dominant correlations 
in this phase has not been clearly worked out. Initally 
Xiang 26 proposed an almost "trimerized" groundstate, 
with dominant singlet correlations involving three con- 
secutive spins. Although models can be constructed 
which have exactly trimerized groudstates^, in the 
present model these are not the dominant fluctuations, 
as shown earlier in Ref. [3(j ■ 

Here we show that the dominant correlations are not 
of singlet, but of spin nematic, i.e. quadrupolar charac- 
ter. We build upon the field theoretical work of Itoi and 
KatcS in which they showed that throughout the critical 
region there are three fields with spin 5 = 0, 1, 2 whose 
scaling dimensions are all equal to x = 2/3. However 
apart from 9 — ir/4 it is the 5 = 2 mode at k = ±27r/3 
- which is subject to correlation- enhancing logarithmic 
corrections - which will show the dominant correlations. 
Note that this is consistent with the observation! that 
the lowest energy level at k = ±2tt/3 on finite chains 
carries 5 = 2. 

In order to work out the correlation content of this 
5 = 2 mode we note that this mode will only show up 
in a dynamical correlation function if it is targeted with 
an 5 = 2 operator. A natural operator is the irreducible 
quadrupolar part of the rank two tensor 5f5f . In order 
to test this conjecture, we calculate a component of the 
static spin quadrupolar structure factor 

Q(k) = ((5 2 ) 2 (-fc)(5*) 2 (fc)) (2) 

and similar for the static spin structure factor 

S(k) = <5 z (~fc)5 2 (fc)) (3) 

both at the wavevector k = 2ir/3. The results shown in 
the upper plot of Fig. for system sizes up to L = 18 
display nicely that the quadrupolar correlations indeed 
become more important than the spin-spin correlations 
when one moves beyond the SU (3) symmetric point at 
9 = 7r/4. Although all correlation exponents (dimer, 
spin, quadrupolar) are equal to r\ — 2x — 4/3, loga- 
rithmic corrections^ render the quadrupolar correlations 
clearly dominant in this phase. Note that our notion 
of a dominant correlation function in this case is analo- 
gous to the 5 = 1/2 antiferromagnetic Heisenberg chain, 
where both staggered spin and staggered dimer correla- 
tions have equal correlation exponent r\ = 1, but the log- 
arithmic corrections enhance the spin-spin correlations 
and then decay more slowly. Hence, in the same spirit as 
one argues that the dominant correlations of the Heisen- 
berg chain are the spin-spin correlations, the dominant 
correlations of the bilinear-biquadratic spin-one chain for 
9 G (+7r/4, vr/2) are of spin quadrupolar nature with 
wavevector ±27r/3. 

For completeness we show selected parameters of this 
critical region in the lower panel of Fig. [21 The results 
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FIG. 2: (Color online) Upper plot: Static structure factors 
for spin and quadrupolar correlations at momentum k — 2n/3 
as a function of 6 for different system sizes. At the SU(3) 
point 8 = +7r/4 the two correlation functions are related by 
symmetry. By going deeper into the critical region the spin 
correlations get weaker, while the quadrupolar correlations 
are enhanced. Lower plot: the central charge, the excitation 
velocity and the scaling dimension of the quintuplet field in 
the critical region. The filled red square is the exact result— 
for the velocity at 9 — 7r/4. 



are in good agreement with field theoretical predictions 2 ^, 
(c = 2 and x q = 2/3), exact results^ [v(B=n/i) = 
27r/3v / 2], and previous numerical work 7 -. The numerical 
estimate of the scaling dimension x q is slightly smaller 
than the expected value of 2/3, due to the presence of 
logarithmic corrections. 

To close this section we note that the occurrence of 
dominant period-three quadrupolar correlations on a sin- 
gle chain calls for a fully ordered three sublattice struc- 
ture when the chains are coupled in an appropriate two- 
dimensional fashion. The triangular lattice perfectly ful- 
fills this prerequisite, and indeed recent worksi 3 ^ re- 
ported analytical and numerical evidence for a three 
sublattice quadrupolar ordered phase in the region 9 € 
(7r/4,7r/2). 



III. FROM THE DIMERIZED TO THE 
FERROMAGNETIC PHASE 

After having elucidated the nature of the dominant 
correlations in the extended critical region beyond the 
Haldane phase, we now turn to the behavior of the 
bilinear-biquadratic chain between the point where the 
dimerized groundstate is firmly established (9 = —w/2)& 
and the point where the ferromagnetic state takes over 
(6 = — 37r/4). Let us note that this point 9 = — 3tt/A is 
special due to a SU(3) symmetry. This symmetry has the 
consequence that the conventional ferromagnetic multi- 
plet is degenerate with a fully ordered ferroquadrupolar 
spin nematic stated. 

Chubukov- suggested the existence of an intermedi- 
ate phase which was supposed to be completely gapped 
and to have dominant short-ranged nematic correlations. 
Since then several studiea 17 i 18 i 19 i 20 i 21 i 22 i 23 - most of them 
numerical - tried to pin down the existence of an inter- 
mediate phase. Presently it seems most likely that the 
gapped nematic phase in its original form is not realized 
in the single chain model, mainly because the gap data 
do not show evidence for a closing and reopening. Never- 
theless the spin nematic correlations grow dramatically 
as one approaches 9 — > — 37r/4 + . It could therefore in 
principle be possible that the dimerized phase gives way 
to a critical spin nematic phase without a gap22i2ii2£. We 
will critically discuss this possibility in the following. 



A. The S — 2 finite size gap 

We track the evolution of the energy gap to the first 
magnetic excitation using Density Matrix Renormaliza- 
tion Group (DMRG) calculations^ on long open chains 
of up to 512 spins and retaining up to 1000 states. 

A phase transition is signaled by the closing of the gap. 
The lowest excited state on open chains carries spin two 
for 9 6 (— 37r/4, — ir/2] 17 . Our results shown in Fig.[31are 
clearly consistent with a finite gap of the 5 — 2 excitation 
for 8 > — 0.657T. However, for 8 < — 0.65tt the extrapo- 
lated gap becomes very small (of the order of 10~ 3 for 
6 = — 0.77r) which could suggest the existence of a phase 
transition around — 0.677T below which the gap is zero. 
The data shows no evidence for a reopening of the 5 = 2 
gap in the interval 8 £ (— 37r/4, — 0.67tt], at variance with 
the initial proposal^. 

In order to investigate the possible closing of the gap 
by different means we have calculated strong coupling 
series expansions 37 of the 5 = 2 single particle gap start- 
ing from the dimerized limit up to 10th order in the in- 
terdimer coupling A for fixed 9. A direct evaluation of 
the quintuplet gap shows a closing of the gap around 
9 = — 0.677T which is illustrated in the left panel of Fig.0] 
Using DLog Pade approximants we have determined the 
critical interdimer coupling X c (9) where the gap closes, 
see upper right panel of Fig. 01 Consistent with previ- 
ous estimates by Chubuko^ these extrapolations suggest 
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FIG. 3: Extrapolated DMRG S = 2 gaps as a function of 
0. The gray region denotes the estimated onset of the critical 
behavior. The filled diamond is the exactly known gap result 
at 6 — — 7r/2. Inset: finite size extrapolation of the S — 2 
gaps (32 to 512 sites). 9/tt varies from —0.5 to —0.7 with 
decrements of 0.025 from top to bottom. 
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FIG. 4: Strong coupling dimer expansion of the S = 2 gap. 
Left panel: S = 2 gap as a function of 8. The gap closes 
at 6 C ~ — 0.67-7T. Upper right panel: Critical value A c where 
the gap vanishes. The uniformly coupled chain of interest 
here corresponds to A = 1 (dashed line). Lower right panel: 
Critical exponent of the S = 2 gap opening as a function of 
the dimerization (A ~ A — A c | 7 ) calculated by various DLog 
Pade approximants. 

that the gap closes for the uniformly coupled (A = 1) 
chain around 9 rs — 0.677T. 

As a further indication for the possible closure of the 
gap around 9 = — 0.677T, we now consider the existence 
of a generalized BKT phase transition. To determine the 
critical point we use phenomenological level spectroscopy. 
We calculate the level crossing 0% of the lowest singlet 
excitation at k — n with the lowest spin-2 level at k = 
for different system sizes L with ED and extrapolate L — > 
oo. The results are shown in Fig.JSJ The extrapolation is 
performed with 1/L and 1/L 2 corrections, fitting only the 



FIG. 5: Finite size scaling of the crossing points between 
the singlet gap at momentum n and the quintuplet (S = 2) 
gap at momentum 0. The extrapolation yields a critical value 
9 C ~ — 0.67-7T in the thermodynamic limit. 



largest four system sizes. The extrapolated critical point 
is 9 C = (— 0.67±0.02)7r, in good agreement with estimates 
obtained by other techniques 2 - 8 . The coefficient of the 
1/L correction is small, and fitting without this term 
yields slightly larger values 9 C > — 0.677T. 



B. Dimerization and quadrupolar correlations 

Dimerization in the phase diagram of the spin-one 
chain has been firmly established by exact results ob- 
tained for the 9 = —ir/2 points. It should be noted that 
although the system is dimerized, it is rather poorly de- 
scribed by a simple product wavefunction of alternating 
singlets, as can be seen from the small gap and the large 
correlation length over the whole dimerized phase. The 
dimerization operator considered here is: 

with k = 7r. The second kind of correlations expected to 
be important in this region of the phase diagramSiiS are 
the ferroquadrupolar spin fluctuations 

Q(A:)^-^^e^[(^) 2 -2/3] (5) 

3 

with k = 0. Note that we considered similar fluctuations 
with k = 2tt/3 in section ILTI 

Based on these fluctuation operators we have calcu- 
lated the standard static structure factor (SF): 

C SF (fc) = ||C(fc)|0)|| 2 (6) 
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FIG. 6: Static and generalized susceptibilities for the bilin- 
ear dimerization (left panel) and the ferroquadrupolar corre- 
lations (right panel), obtained by ED on systems of 8 to 16 
sites. Solid symbols: Static structure factors. Open symbols: 
Corresponding generalized susceptibilities. 



and the generalized nonlinear susceptibility (i.e. the real 
part of the dynamical correlation function at zero fre- 
quency) 



tGNS 



(fc) = H lim MC\-k)— 



IT] 



-C(k)\0) (7) 



for both kinds of correlations (C = T>, Q). The nonlinear 
susceptibility quantifies the perturbative response of the 
system upon explicitly coupling the symmetry breaking 
operator to the Hamiltoniari2&. 

We present ED results concerning the dimerization in 
the left panel of Fig.© The SF and the GNS both diverge 
as L — > 00 for 9 = —tt/2, where long range dimer order is 
established 6 . Similar behavior is found for 9 > — 0.677T. 
For 9 e (-37I-/4, -0.67tt) however, we find that the SF 
and the GNS both decrease with system sizes, pointing 
to a possible absence of spontaneous dimerization in this 
region. 

In the right panel we display the same kind of ob- 
servables for the k — (ferroquadrupolar) mode of the 
quadrupolar correlations. Here the behavior seems dif- 
ferent: while there are only short range correlations deep 
in the dimerized phase, the ferroquadrupolar correla- 
tions increase drastically with system sizes for 9 close 
to -3tt/4. 

In order to shed more light on the excitations and their 
nature close to 9 = —3tt/A we have calculated in addition 
the dynamical dimer and quadrupolar structure factor 
using a continued fraction technique. This will allow us 
to track the evolution of the energy of the important 
levels and their spectral weight. The dynamical structure 
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FIG. 7: ED dynamical correlation results on a L — 16 chain. 
Upper panel: dynamical staggered dimer correlation function 
T>(q — 7T, u)) plotted for a range of values of 6. There is an 
important transfer of spectral weight upon lowering 8 from a 
low energy level - collapsing to the groundstate as a function 
of system size - to a level which will converge to the energy 
\/2 at 6 = — 37r/4. Lower panel: dynamical ferroquadrupolar 
correlations Q(q = 0,u>) for 9 G (0.75, 0.60]tt. 



factor is defined as follows: 



C(k,u) 



= lim (0\CU-k) 

7T ri~*0+ 



1 



(H-Eo 



- ir] 



-C(fc)|0>, 

where we focus again on the dimerization and the fer- 
roquadrupolar correlations. The plots discussed below 
show the normalized intensity for a given 9. The overall 
weight can be obtained by multiplying with the structure 
factors displayed in Fig. 

The results for the dimerization shown in the upper 
panel of Fig. [7| are very interesting. They show a strong 
crossover with a redistribution of spectral weight on the 
way from 9 = —tt/2 to — 3tt/4. While there is a single 
state (the lowest energy singlet) exhausting almost all of 
the k = it dimer weight around 9 = —tt/2 (as expected 
in a truly dimer ordered phase), the weight in this state 
entirely fades out and is transfered to a higher energy 
state at uj — > %/2 as 9 — > — 37r/4. 
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FIG. 8: (Color online) Hypothetical conformal field theory 
parameters calculated in ED on chains of 8 up to 18 sites. The 
scaling dimensions of the singlet (k — n), and the quintuplet 
(k = 0) field are shown. Inset: effective central charge. The 
value at the boundary of the critical region (8 « — 0.677r) is 
cw3/2. 



The results for the dynamical ferroquadrupolar corre- 
lations don't show such a crossover but show a steady- 
lowering of the finite size S = 2 gap upon approach- 
ing — 37r/4, accompanied by an accumulation of all the 
weight in that single state. This is a manifestation of the 
continuity of the finite size groundstate wavefunction as 
one approaches the fully ferroquadrupolar ordered state 
at 9 = — 37r/4. We return to the special behavior of the 
spin- two gap in subsection 1111 C 2l 

We believe that it is this strong crossover seen in the 
dynamical dimer correlations which renders numerical 
calculations very difficult. On finite size systems suffi- 
ciently close to 9 = — 3ir/4 the system looks basically 
like a ferroquadrupolar ordered system, even in one di- 
mension, due to the very strong influence of the SU(3) 
point — 37r/4. As a consequence a crossover scale seems 
to emerge, which rapidly grows close to — 37r/4, and one 
has to go to huge systems in order for the dimerization 
to win, i.e. for the lowest singlet at k — ir to collapse 
onto the groundstate, while the lowest 5 = 2 exitation 
at k — should remain gapped. This scenario would 
then imply that the curve in Fig. [5] finally would have to 
bend down for large systems and touch the y-axis only 
at 9 = -3tt/4. 



C. Critical phase versus Crossover 

1 . Parameters of a hypothetical critical theory 

In order to discriminate between potential conformal 
theories describing the phase transition and the extended 
gapless region, we calculate relevant field theory param- 
eters, i.e. the scaling dimensions (SD) x of several fields 
and the central charge c (Fig. [HJ) - The calculation of these 
parameters has been performed with ED and relies on fi- 



nite size scaling properties of groundstate and excited 
states energies^. Leading logarithmic corrections have 
been taken into account. The SD of the k = ir singlet 
field (x s ) and the k = quintuplet field (x q ) are both 
close to 3/8 at the onset of criticality near 9 w — 0.677T. 
x s significantly increases, while x q decreases and seems 
to approach as 9 — > — 37r/4, compatible with quadrupo- 
lar long range order precisely at 9 = —3ir/4£^. Our re- 
sults for the effective central charge are unexpected. The 
central charge c does not seem to be constant through- 
out the potentially critical region. The smallest value is 
found at the onset (c ~ 3/2), and then seems to increase 
monotonously. While this behavior can not directly be 
ruled out on field theoretical grounds, it is rather un- 
common. It remains to be understood whether c is re- 
ally continuously increasing or whether we are facing a 
crossover phenomenon. The hypothetical critical theory 
at 9 rs — 0.677T is however surprisingly well characterized 
by a level two SU(2) Wess-Zumino-Witten model given 
that both the scaling dimensions (3/8) and the central 
charge (3/2) are in agreement with such a theory. Some 
more support for this claim comes from the indirect cal- 
culation of x s within our series expansions. There the 
critical exponent 7 of the S — 2 gap is related to the 
scaling dimension a: s by 7 = 1/(2 — x s ). In the lower 
right panel of Fig. 0] we show the critical exponent cal- 
culated by DLog Pade approximants to the 10th order 
series. In the vicinity of 9 rs — 0.677T the various approxi- 
mants reveal only small spreading. Within the precision 
of a 10th order calculation the exponents comply with a 
scaling dimension of x s =3/8 (dashed line). 



2. Single Mode approximation 

In the following we present an argument based on the 
single mode approximation (SMA) which aims at explain- 
ing the anomalously strong suppression of the S = 2 gap 
upon approaching — 37t/4. 

In the single mode approximation one starts by con- 
structing a trial state upon the application of a structure 
factor operator on the groundstate: 



|0(*)> = 0(A)|O> , 



where 



0{k) 



Akr 



id 



(9) 



(10) 



and Oj is an operator acting on site j. One can now relate 
the variational energy of the state \0(k)) to a groundstate 
expectation value: 



l {0\[OH-k),[H,O(k)}}\0) 
2 (O|0t(-*)0(fc)|O) 



(11) 



Note that the energy lo^ is a strict upper bound on the 
gap in the momentum k sector. The power of the SMA 
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comes from the fact that it can be used to prove the ab- 
sence of a gap under certain conditions. There are two 
prototypical cases for the absence of a gap: i) The de- 
nominator diverges as a function of systems size, while 
the numerator diverges more slowly or stays finite. This 
is the conventional situation for systems which are crit- 
ical or exhibit spontaneous symmetry breaking in the 
thermodynamic limit. In this case the diverging static 
structure factors drives the gap to zero, but only for the 
infinite system, ii) The numerator vanishes, while the 
denominator does not vanish. A sufficient condition for 
the numerator to vanish is if the structure factor operator 
0(k) commutes with the Hamiltonian. 

If we now choose for the operator O(k) the spin 
quadrupolar structure factor at zero momentum: i.e. 
Q(0) = l/VLJ2i [(S z )t - 2/3] we realize scenario ii) at 
the special SU(3) point 8 = — 3ir/4. For this specific 
value 2(0) commutes with the Hamiltonian and there- 
fore the numerator is zero. Numerically we find that the 
denominator is different from zero. This means that there 
is at least a level with 5 = 2 which is degenerate with the 
groundstate already on finite size samples^. The reason 
which makes the numerical analysis for 8 > — 37r/4 so 
difficult is that there is a continuity in the wavefunction 
as a function of 8 when approaching the 5t/(3) point 
from above. On one hand we know that the commuta- 
tor [2(0), H(8)] will vanish continuously as we approach 
— 37r/4 + , on the other hand we find numerically that 
the denominator (the ferroquadrupolar structure factor) 
grows rapidly as a function of the accessible system sizes, 
c.f. full symbols in the right panel of Fig[5] So these two 
factors cooperate in giving an anomalously small SMA 
energy u>k, which is an upper bound to the measured 
gap. 

IV. TWO COUPLED CHAINS 

We have seen considerable difficulty in the preceding 
section to actually decide whether the dimerization van- 
ishes before reaching 8 = — 37r/4. One potential way 
of circumventing the spontaneous dimerization of a sin- 
gle chain is to consider two coupled chains in a ladder 
geometry. It seems rather natural to assume that the 
dimerization of a single chain disappears once the cou- 
pling on the rung is suffiently strong. We will substan- 
tiate this claim using numerical simulations below. The 
ladder model is also a first step towards a realistic setup 
of spin-1 bosonic atoms in optical lattices, where a finite 
interchain coupling can easily be generated. The ladder 
Hamiltonian is given as follows: 

H = J ^ cos 8 (Si, n ■ Sj+i, w ) + sin 8 (Sj |Tt • S i+1 , n ) 2 

i,n 

+K cos 9 ( S -i • + sin 8 (Si,! • S 4 , 2 ) 2 , (12) 

i 

where Sj in denotes a spin-one operator at position i on 
chain n G {1,2}. In the following we choose K, J > and 
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FIG. 9: (Color online) Phase diagram of two coupled bilinear- 
biquadratic 5 = 1 chains. The dimerization along the chains 
is rapidly suppressed by a finite interchain coupling. The 
dominating phase is a unique gapped phase which is adiabat- 
ically connected to the Haldane phase of two isolated Heisen- 
berg chains. 



vary 8 as well as Kj J > 0, where the limit K/J = cor- 
responds to decoupled chains and J/K = to decoupled 
rung dimers. 

For 8 = the Hamiltonian l|12fl describes a conven- 
tional 5=1 Heisenberg spin ladder 25 40 . Interestingly, 
for this special case Todo et alm^ have shown numer- 
ically that the finite gap of the single chain does not 
close for any K/J, i.e two points on the Kj J axis can be 
reached without any second order quantum phase tran- 
sition in between. They also generalized the string order 
parameter of a single chain to a more complicated non- 
local operator which is non-zero for any finite value of 
K/J, therefore showing that this phase is also topologi- 
cally ordered 25 . In the following we will show that this 
gapped spin liquid state at 8 = extends deeply into 
the region 8 < 0, and for suffiently large K/J even to 
— 37r/4 + . While the groundstate is always protected by 
a finite gap, the nature of the lowest excitation changes 
as a function of 8. For 8 € [— tv/2,0] the lowest excita- 
tion carries 5 = 1, while for 8 6 [— 37r/4, — tt/2] it carries 
5 = 2. In the absence of dimerization such a state on the 
ladder where the lowest excitation is of ferroquadrupolar 
nature, is very close in spirit to the long-sought gapped 
nematic phase initially proposed by Chubuko"^ for the 
single chain. The phase diagram of two coupled chains 
with 8 e [-7T, 0] and K/J £ [0, 1] is shown in Fig. 
We expect the extended gapped phase to continue up to 
the rung dimer limit J/K = 0, similar to the results for 

= 025,40. 

Let us now describe the simulation results which lead 
to the phase diagram presented in Fig-EI First we discuss 
the behavior of the dimerization of the single chain upon 
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FIG. 10: (Color online) Finite size scaling of the dimerization 
measured in the middle of a ladder as a function of K/J. It 
can be clearly seen that that the dimerization extrapolates 
to a finite value for K/J < 0.125, while it tends to zero for 
K/J > 0.15. 



FIG. 11: (Color online) Evolution of the spin gap for constant 
K/J = 1 as a function of 6 £ [— Sn/4, 0]. The nature of the 
gap changes at the SU(3) point 8 — —it/2, from a S=l state 
for 8 > -it/2 to a a S=2 state for 8 < -tt/2. The QMC value 
of the gap at 8 — is from Ref. |25|. 



coupling two chains. To this end we performed DMRG 
calculations on ladder systems of sizes up to 2 x 256 sites. 
We have determined the boundary of the dimerized phase 
using finite size extrapolations of the remnant dimeriza- 
tion in the middle of a chain^i, assuming a columnar 
dimer arrangement on the two chains^. An example for 
such finite size data is shown for 9 = —tt/2 in Fig. ITUl 
We expect the dimerization in the middle of the chain 
to converge to a finite value for 1/L — > in a dimerized 
phase, and to exhibit an exponential drop upon reaching 
the correlation length of a disordered phase. These two 
distinct behaviors can be seen in Fig. ^| for the values 
K/J = 0,0.1,0,125 and K/J = 0.15,0.175,0.2,0.3 re- 
spectively. This leads us to the conclusion of a critical 
ratio {K/J) c = 0.1375±0.0125 at = -tt/2, below which 
the ladder is still dimerized. The evolution of the values 
of (K/J) c with 9 is shown in Fig.0 The dimerized phase 
is most stable close to 9 — —tt/2, and seems to vanish 
linearly upon approaching the BT point at — tt/4, similar 
to the case of explicit dimerization for a single chainii 
On the other side however the dimerization boundary 
drops rapidly and becomes very small once 9 < — 0.77T. 
Similar to the single chain case it is very difficult to de- 
cide whether the dimerized phase ceases to exist before 
— 37r/4 or not. In the ladder case however this implies 
that for a finite - but very small - interchain coupling 
the spontaneous dimerization (if any) of a single chain 
vanishes. 

Next we investigate the behavior of the spin gap in the 
non-dimerized phase. We performed DMRG calculations 
for a fixed value of K/J = 1. The results for the gap to 
the lowest S z = 1 state are plotted in Fig. ^2 for system 
sizes 2x1 with L = 32, 64, 128. Starting at 9 = where 



we find good agreement with the high-precision QMC 
gap from Ref. |25j, the gap grows monotonously until 
reaching the bipartite SU(3) point at 9 = —tt/2. There 
we observe a level crossing, and the gap decreases for 
9 < —tt/2. Note that at this special point the total spin 
quantum number of the lowest level changes. While it is a 
5 = 1 state for 9 > —tt/2, the lowest levels carries 5 = 2 
for 9 < —tt/2. Similar to the single chain case discussed 
in Sec. IIII Al we detect a very small gap for 9 close to 
— 37r/4. This anomalously small 5 = 2 gap can again be 
understood using a single mode approximation argument, 
along the lines of Sec. 1111 C 21 the extension of the argu- 
ment to the ladder case being straightforward. We take 
this dichotomy of the gap as a first indication that the 
SU(3) line (9 = —tt/2) constitutes a crossover line inside 
the extended gapped phase, where the spin-spin correla- 
tions dominate for 9 > —tt/2, while we expect dominant 
short-ranged spin nematic ferroquadrupolar correlations 
for 9 < —tt/2. Our gap results give evidence that the 
groundstate for 9 close to — 37r/4 adiabatically connects 
to the "Plaquette Singlet Solid" phase, which is itself 
connected to the Haldane phase of the isolated chains. 
This would also imply that the the whole phase is char- 
acterized by a finite nonlocal order parameter. It would 
be interesting to investigate this order in a future work. 



V. CONCLUSIONS 

To summarize, we have shown that spin nematic cor- 
relations in the form of quadrupolar correlations play an 
important role in the phase diagram of single or coupled 
5 = 1 bilinear-biquadratic chains. 
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We first addressed the long-standing question of the 
nature of the dominant correlations in the gapless period- 
three phase, where we uncovered the spin quadrupolar 
correlations at k = ±2tt/3 as the leading ones. We crit- 
ically discussed the phase diagram close to the SU(3) 
point at 6 = —in/ A and by assembling several numerical 
and analytical results concluded that an unconventional 
crossover phenomenon is at the heart of the consider- 
able difficulty in settling the issue on the existence or 
absence of the gapped nematic phase put forward by 
Chubukov 8 . Finally we studied two coupled bilinear- 
biquadratic chains in a ladder geometry and found an 
extended gapped phase which is surprising in two re- 
spects: first it realizes a variant of Chubukov's nematic 
phase, i.e. a gapped, non-dimerized phase with dominant 
spin-nematic correlations for < —n/2, and second, this 
phase is adiabatically connected to the standard 5 = 1 



Heisenberg (9 — 0) ladder, and therefore also connected 
to the Haldane phase of single chains, as shown by Todo 
et al&. 
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